Submitted to ApJL 



Expansion Velocities and Core Masses of Bright Planetary 
Nebulae in the Virgo Cluster 

Magda Arnaboldi^'^, Michelle Doherty^, Ortwin Gerhard^, Robin Ciardullo^ 

00 " 

^ ! J. Alfonso L. Aguerri^, John J. Feldmeier^, Kenneth C. Freeman'', George H. Jacoby^ 

g : ABSTRACT 



The line-of-sight velocities and [OIII] 5007 A expansion velocities are mea- 
sured for 11 planetary nebulae (PNs) in the Virgo cluster core, at 15 Mpc dis- 
Qh' tance, with the FLAMES spectrograph on the ESO VLT. These PNs are located 

O ■ about halfway between the two giant ellipticals M87 and M86. From the [OIII] 

^ • 5007 A line profile widths, the average half-width at half maximum expansion 

velocity for this sample of 11 PNs is vhwhm = 16.5 kms~^ (RMS = 2.6 kms~^). 
We use the PN subsample bound to M87 to remove the distance uncertainties, 
^ ■ and the resulting [OIII] 5007 A luminosities to derive the central star masses. 

^ ■ We find these masses to be at least 0.6 M© and obtain PN observable life times 

O ■ tpiy < 2000 yrs, which imply that the bright PNs detected in the Virgo cluster 

core are compact, high density nebulae. We finally discuss several scenarios for 
^ ■ explaining the high central star masses in these bright M87 halo PNs. 
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Introduction 



Sin ce the discovery of free -floating; intracluster planetary nebulae (ICPNs) in the Virgo 
cluster ( lArnaboldi et al.lll996l ). extensive imaging and spectroscopic observations were car- 
ried out to determine their projected phase space distribution and the fraction of diffuse 
cluster light not bound to Virgo galaxies. To enlarge the s ample of rn ore than 40 ICPN 



line-of-sight (LOS) velocities available from lArnaboldi et al.l (2003,1200411 . we have obtained 



new PN spectra with FLAMES at the ESO VLT flDohertv et al. 



20081 ). The high spectral 



resolution of the new data allows us for the flrst time to measure expansion velocities for 
the [OIII] nebula shells of 11 PNs in the Virgo cluster core, nearly half way between M87 
and M86. The expansion velocity of a planetary nebula (PN) is one of the most impor- 
tant parameters determining its evolution, but currently it is known only for a few' hundred 
Galactic PNs, mostly bright objects in the Milky Way Disk (iGesicki &: Zijlstrall2000l ). and for 
a few tens in the Ma gellanic Clouds and M31 . When interpreted using dynamically evolving 



nebular models (e.g. Schonberner et aTlboOsI ). PN shell expansion velocities provide reliable 



extimates of PN dynamical ages and distance estimates for Galactic PNs. 

In this letter, we give a brief summary of our Observations in Section [2l In Section [3l 
we present the PN [OIII] 5007 A half-w idth ha l f max imum velocity measurements vhwhm-, 
and m(5007) magnitudes as deflned by I Jacobyl (119891 ). In Section H] we build the PN lumi- 
nosity function (PNLF) for the spectroscopically conflrmed PN sample in Virgo and for the 
subsample bound to the M87 halo. In Section O we then estimate the outer radii of the 
M87 halo PNs from the observed PN expansion velocities and their visibility lifetimes. We 
flnally discuss in Section [6] the remaining distance uncertainties, the shape of the PNLF in 
the M87 halo, and the mechanisms that may be responsible for the large core masses of the 
brightest PNs in the halo of M87. 



Observations 



Data were acquired in service mode (22 hrs, 076.B-0086 PI: M. Arnaboldi), with the 
FLAMES spectrograph at UT2 on VLT in the GIRAFFE+MEDUSA conflguration. We 
used the high-resolution grism HR 504.8, covering a wavelength range of 250 A centered 
on 5036 A and a spectral resolution of ~ 20000. The redshifted [OIII] emissions of PNs 
in the Virgo cluster core fall near the center of the grism transmission. With this setup, 
the instrumental broadening of the arc lines is FWHM = 0.2 9 A or 17 kms~\ and the 
error on the wavelength measurements is 0.0025 A or 150 ms~^ ( Royer et al.ll2002l ). A total 
of three FL AMES plate con fl gurat ions were produced for the Virgo core flelds F4 and F7 
surveyed by iFeldmeier et al.l (120031 ). The exposure times were based on the signal-to-noise 
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ratio (S/N) estimate for detecting the [OIII] 5007 A line flux of 4.2 x 10"^''' ergs cm~^ s~\ 
i.e., m(5007) = 27.2, with S/N ~ 5. The data reduction was carried out with the GIRAFFE 
pipeline, for the CCD prereduction, fiber identification, wavelength calibration, geometric 
distortion corrections, co-addition, and extraction of the final one- dimensional spectra. Data 
analysis and the velocity measurements are presented in Doherty et al. (2008). A total of 
12 PNs were confirmed spectroscopically. In the extracted PN spectra, the [OIII] 5007 A 
emissions of 11 PNs have measured FWHM between 0.4 and 0.7 A; the comparison with 
the archnes' FWHM (0.29 A) indicates that these [OIII] emission lines are resolved. Monte 
Carlo simulati ons give an error on the line width measurements of 0.0025A for a S/N per 
pixel of ~ 15 jRover et al.ll2002h . 



3. PN expansion velocities 



The spectroscopic expansion velocity measured from the [OIII] 5007A line profile ap- 
pears to be representative for the material velocities associated with the PN bright ri and 
is systematically lower than the ex pansion speed of the PN shell's outer radius, according 
to ID hydrodynamical simulations (ISchonberner et al.ll2005l ). The shell expansion velocities 
are associated with fainter struct ures at larger radii, which have been measured only for 
Galactic PNs dCorradi et aDl2007h . 



From the FLAMES PN spectra, we are able to measure the sp ectroscopic expansion 
velocities from the resolved [OIII] 5007A emissions for 11 PNs in the iDohertv et all (120081 ) 



sample; high S/N spectra for 4 PNs are shown in Fig. [H We use the "half- width at half 
maximum" of the [OIII] 5007 A line, corrected for the instrumental half width, vhwhmi as 
measurement of th e spec trosc opic expansion veloci ty for the PN bright rim. According to 
Schonberner et al. J2005h and ICorradi et al I (120071 ). these vhwhm measurements must be 
multiphed by a factor of about two to get an estimate of the PN's true expansion velocity. 
When it is possible, we also measure the spectroscopic expansion velocities for [OIII] 4959 A. 
When double peaks or secondary peaks are present in the resolved line profiles, their counts 
are consistent with the noise. The average [OIII] spectroscopic expansion velocity for this 
PN sample is vhwhm = 16.5 kms~^, with an RMS dispersion of 2.6 kms~^. Both values 
are sig nificantly smaller than t hose determined for samples of PNs observed in the Galactic 
bulge (iGesicki fc Zijlstrall2000l ). The Virgo PNs might represent the [OIII] brightest part 
of the Galactic bulge population, or they might come from an entirely different population; 
this issue requires further work. 



The bright rim is the thin shcU weU behind the outer shock, enclosing the wind-blown cavity of a PN. 
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Based on LOS velocities, we divide the current dataset in two subsamples: 5 PNs 
are associated with the halo of M87, i.e. belong to a narrow velocity peak in the LOS 
velocity distribution (LOSV D) centered at the M87 systemic velocity, with RMS = 78 



kms ^ (jPoherty et al. 



20081 ) ■ and 6 are "free-flying" cluster PNs, with 300 - 1500 kms' 



velocity differences. The distributions of the vhwhm measurements for the two subsamples 
are similar within the limited statistics, see Fig. [21 the KS-test gives only a 55% probability 
that the two subsamples are different. The vhwhm vs. m(5007) plot in the top panel of 
Fig. [2] shows that the PNs within 0.7 mags of the bright cut-off have vhwhm < 20 kms^^. 
This is in agreement with the predictions of the spatially integrat ed line proflle vhwwm vs - 
m(5007) computed with the nebular ID hydrodynamics code of ISchonberner et al.l (120071 . 
and in prep.) for central star masses in the range 0.696 — 0.625 Mq. 

When a PN's distance is known, its expansion velocity and size can be used to determine 
its dynamical age. This age, in turn, can be combined with infor mation about the centra l 
star's effective temperature to yield an estimate of the core mass ( iGesicki fc Zijlstral 120071 ). 
In our case we cannot measure the outer radii for the Virgo PNs, as they are unresolved at 
a distance of 15 Mpc. Our approach is then the reverse: we shall estimate the central star 
masses from the [OIII] 5007 A fluxes, derive the PN lifetimes tpj^f from central star evolu- 
tionary tracks and nebular 1-D hydrodynamical models, and infer their physical dimensions 



as rpAT 



X tpisf, where Vexp = 2 x vhwhm (Schoenberner et al., in prep.) 



4. The M87 halo PNLF 

The PNLF for the spectroscopically conflrmed PNs in the Virgo core region around 
M87 is shown in Figure O The brightest PNs have m(5007) in the range 25.7 — 26.5. PNs 



with similar magnitudes were also de tected bv ICiardu 



lo et al.l ( Il998l ) in the outer regions of 



M87. Based on a sample of 338 PNs, ICiarduUo et al.l (119981 ) determined a brightening of the 
PNLF cutoff of 0.37 mag for the PN subsample at i? > 4' with respect to the sample inside 
R < 4'. They interpreted the brightening of the PNLF in the M87 outer halo as due to a 
population of Virgo ICPNs flUing the elongated volume of the Virgo cluster (up to 4 Mpc 
along the LOS). The number of foreground PNs would be roughly proportional to the area 
of the fleld, and therefore be largest in the outer regions of the surveyed fleld. 

In our spectroscopic PN sample we can bypass the distance ambiguity by selecting the 
subsample of 14 PN^O bound to M87. These PNs are associated with a narrow peak in 
the LOSVD centred at the systemic velocity of M87, see Fig. |3l Their average velocity 



2 5 PNs from the lOohertv et al.l |2008l) sample and 9 from the F3 field from lArnaboldi et all |2004h 
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vlos = 1306 kms~^ and the RMS = 117 kms~^. The empirical PNLF for the M87 halo 
PNs is shown in Fig. [3j their m(5007) magnitudes are in the range 26.2 — 27.2, indicating 
a slightly brighter cut-off than for the PN population in the central R < A' region of M87. 
However we can now say that this is an intrinsic property of the PN halo population, as we 
now know that these 14 PNs are all at the distance of M87 within ~ 100 kpc (see also §6). 



5. PN core masses, visibility time scales and outer radius 



From the measured m(5007) of the M87 halo PNs, we estimate the central star total 
luminosity. Both models and observations indicate that no more than 10% of the central 
star's total luminosity comes out in this line (jJacobylll989l ). Therefore the intrinsic luminosity 
of the post-AGB star that powers an m [5007) = 26.2 P N at 14.5 Mpc must be > 6930 L© 
and, based on the evolutionary tracks of iBlocker I (119951 ). we thus obtain a central core mass 
l arger than 0.6 M^)- Using a 1-D hydrod ynamics code and the stellar evolution tracks of 
( iBlocker I Il995l ) , ISchonberner et al.l (120071 ) computed the nebular physical parameters and 
evolution of m(5007) as function of nebular age, from near th e AGB phase to the white 
dwarf cooling tracks. The nebular tracks of ISchonberner et al.l (120071 . Fig. 15) that reach 
the brighest 0.5 mag of the PNLF, and have a central core mass of > O.6OM0 as the Virgo 
PNs, have very short visibility lifetimes, tp^ < 2.0 x 10^ yrs. From our measurements and 
Vexp = 2 X vhwhm, we can then determine the outer radii rpjy of the brightest PNs in the 
M87 halo to be ~ 0.07 pc. The nebular shells of these PNs are compact and similar to those 
observed for the brightest (logL > 3.8) Galactic Bulge PNs. 



Discussion 



Distance uncertainties and PNLF brightening - We briefly consider the question whether 
the bright PNe in the velocity range bound to M87 could be foreground obje cts. The mean 
heliocentric velocity of the Virgo cluster is < Vq >vc= 1050 ± 35 kms~^ (IBinggeli et al. 
19931 ). M87 is redshifted by about 300 kms~^ {vm87 = 130 7 kms~^) with respect to the 



Virgo cluster mean velocity, and is falling into Virgo from in front (IBinggeli et al.l Il993l ) 
towards M86 JPohertv et al.ll2008[ l. 



PNs at the cutoff m* = 26.33 with apparent magnitudes 26.2-26.0 would be between 
~ 1 — 3 Mpc in front of M87. Diffuse light stars at these locations would be in the infall 
region towards Virgo where the infall velocities are of order 1000 kms~^ and vary rapidly 
with distance (iMohayee fc TuUy I l2005l . and references therein). On the other hand, radial 
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velocities of ~ 1300 kms~^ would be expected for foreground objects ~ 6 Mpc in front of the 
Virgo core. It is thus very unlikely that a distribution of foreground PNe should be observed 
in our fields at exactly the systemic velocity of M87, with a dispersion of only ~ 100 kms~^, 
and less likely still that this population should be projected onto the M87 halo inside 150 
kpc but not be observed in the adjacent Virgo core region covered by our data. 

While the Virgo cluster has a significant depth (~ ±2Mpc), it is unlikely that the 
diffuse light and ICPN distribution are similarly elongated. The search for ICPNs outside 
the Virgo clust er core has g i ven n egative results (Castro et al. 2008, in prep.), and the deep 
photometry of iMihos et al.l (120051 ) shows that the diffuse light is mostly associated with the 
giant elliptical galaxies, M87, M86 and M84, whereas its surface brightness decreases sharply 
at larger radii. Observations of the diffuse light in z = 0.2 5 galaxy clusters al so show that it 
is more centrally concentrated than the cluster g alaxies (IZibetti et al.ll2005l ). in agreement 
with cosmological simulation of cluster formation ( Murante et al. 20071 ). Thus we can safely 
conclude that the selected M87 halo PNs are bound to M87, and that their bright [OIII] 
magnitudes are an intrinsinc property of this PN population. 

Mechanism leading to large core masses in the M87 halo PNs - The initial mass-final 
mass relation (IFMR) for solar metallicit y stars predicts that the PN progenitors w ith ~ 2.2 



CiarduUo et al. 2005 



Buzzoni et al. 12006 ). Turnoff 



Mq give final core masses of 0.62 Mq (see 
masses of ~ 2 Mq belong to populations with ages ~ 1 Gyr (llben &: Laughlinlll989l ). The 
question is whether such populations exist in the M87 halo or in the Virgo diffuse stellar 
light. Observations in the V, I bands of the stellar population in a 11.39 arcmin^ field halfway 
between M87 and M86 were carried out with the Ad vanced Camera for Sur vey (ACS) and the 
HST. This location f alls within the F4 Virg o field of lFeldmeier et al.l (120031 ) and is included in 



the area surveyed by lDoherty et al.l (120081 ). At this location. IWilliams et al.l (120071 ) detected 
some ~ 5300 intracluster red giant branch stars (IRGB); from the color magnitude diagram 
(CMD), they estimated the age and metalhcity distribution of the parent stellar population. 
In this region, 70% - 80% of the Virgo IRGB are old (> 10 Gyr), and span a wide range of 
metalhcities (-2.3 < [M/H] < 0.0), with a mean value of [M/H] ~ -1.0. 

From the number of PNs within 0.5 mag of the PNLF cut-off, we can determine 
how much luminosity would be present in an intermed i ate ag e population and compare it 
wit h the measurec 
the Williams et al 



surfac e photometry of iMihos et al.l (120051 ) in the surveyed region and 
([20071 ) fit to the IRGB CMD. If we assume the analytical formula of 



CiarduUo et all (Il989t l for the PNLF, then iVo.5 = Npn/100, where Np^ is the total number 
of PNs associated with the luminosity of the parent stellar popul ation. The luminosity- 
specific PN number a = Np^/Lgai is given by iBuzzoni et al.l (120061 ) for stellar populations 



with different ages and metalhcities, and calibrated using the PN population in the Local 
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Group, Leo group, and the Virgo and Fornax clusters. The maximal theoretical value of 
a which gives the largest PN population for a given luminosity is oimax = ^PN x (1.85 x 
IO^Lq)"^. This theoretical maximal value is independent of metallicity, and a lso provides an 



uppe r limit to the observed a for the PN population in different galaxy types (IBuzzoni et al. 

mm- 



In the M87 halo, we have Ni 



0.5 



minimal bolometric luminosity of the parent population of Lgai = Npjq/amax = 1-66 x IO^Lq 



9 ±3, which gives a total population of A^pat = 900 and a 
en 

over a total surveyed area of 570 arcmin^. We can thus derive a lower limit to the mean 
surface brightness in the V band, /ly = 28.6 mag arcsec"^, for a possible intermed iate age 
parent popula tion. Comparing wit h /xy = 28.3 measured by iMihos et al.l (120051 ) at the 
position of the lWilliams et al.l (120071 ) field, we find that, to justify the number of bright PNs, 
at least 50% of the stars in this field would need to come from a 1 Gyr population. However, 
the uppe r limit to the contribu tion of a younger (< 10 Gyr) component to the IRGB stars 
given by I Williams et al.l (120071 ) is 3 0% - 20% . We must then conclude that there are not 
enough intermediate age stars in the lWilliams et al.l (120071 ) field to justify the observed N^,^. 



What are the pos sible alternatives, if an intermediate age population is not present? The 
Williams et al.l (120071 ) results indicate that the Virgo core population of stars is dominated 
by low metallicity stars ([M/H]< —1) with ages > 10 Gyr; thus we may argue that the halo 
of M87 may be also metal poor. In this case, the PNLF might have a brigh t cutoff that 
is different from that of a metal rich population. However iDopita et al.l (119921 ). and PNLF 
observations in metal-poor galaxies (jCiarduUo et al.l |2002|) show that metal poor systems 
have values of M* that are fainter than those of t heir metal-rich counterparts. Also the 
Oxygen abundance measurements by lMendez et al.l (120051 ) of the brightest extragalactic PNs 
in NGC 4697 indicate near solar metallicities for these stars. Therefore PN evolution from 
a metal poor population is unlikely to be a viable explanation for the observed large core 
masses in the M87 halo PN. The evolution of single stars from a 10 Gyr old p opulation leads 
to central star masses in the range 0.52 <M< 0.55 (IBuzzoni et al.ll2006l ). which cannot 
sup ply the fOIII] 5007A flux required at the PNLF bright cut-off of M* = -4.48. This 
led ICiarduUo et al.l (120051 ) to propose an alternative form of evolution, i.e close binaries and 
blue stragglers stars, as the likely progenitors of [OIII] - bright PNs in a 10 Gyr old stellar 
population. This evolutionary channel seems better in agreement with the observations for 
the brightest PNs in the M87 halo, than either the ~ 1 Gyr-old or the metal-poor progenitors. 
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7. Conclusions 

We have measured the nebular [0111] 5007A spectroscopic expansion velocities for 11 
PNs in the Virgo cluster core, at 15 Mpc distance, with the FLAMES spectrograph on 
the ESO VLT. Based on the [0111] line profile width, the average spectroscopic expansion 
velocity for this sample is vhwhm = 16.5 kms^^ (RMS = 2.6 kms~^), which is in agreement 
with the predictions of dynamically evolving nebular models for high density nebulae close 
to their maximal m(5007) emission. Large central star masses Mcs > 0.6 Mq are inferred 
from the bright measured [0111] luminosities and the known distance for the PNs bound 
to the M87 halo. From the large central star masses and the measured vhwhm, "we derive 
short PN visibility times and small nebular outer radii, ~ 0.07 pc. The PNs in the M87 halo 
have large central star masses, are compact and their nebula shells may be similar to those 
observed for the brightest Galactic Bulge PNs. Three mechanisms are reviewed as possible 
explanation for the large core masses of the M87 halo PNs: intermediate age population, 
metallicity effects, and blue stragglers, with the latter being the most likely, given the old 



age and the low metallicities of the IRGB stars in the Virgo core (jWilliams et al.l 120071 ) 



We thank the referee, Dr Detlef Schonberner, for sharing his results before publication 
and for his constructive comments. We thank Nando Patat and Marina Rejkuba for support. 
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Fig. 1.— Resolved [OIII] 5007 A emission lines of 4 spectroscopically confirmed PNs. 
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Fig. 2. — Top panel - [OIII] spectroscopic expans ion velocities vhwhm vs. m(5007) magni 
tudes for the PNs from the iDoherty et al.l (120081 ) sample; the circled symbo ls indicate the 
M87 P Ns. Bottom panel - histograms of vhwhm measured for all PNs in the iDoherty et al. 



(120081 ) sample (dashed blue hues), for the subsample bound to the M87 halo (continuous red 
lines), and for the PNs not bound to M87 (dash-dotted green lines). The average vhwhm 
for the entire sample is 16.5 kms^^. 
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Fig. 3. — Left panel: PNLF in the surveyed fields in 0.3 mag bins. Long dashed blue line: 
PNLF for the entire spectroscopically confirr ned sample of PNs in the Virgo cluster core 
at the position of the F3, F7 and F4 fields of iFeldmeier et al.l (120031 ) . Red continuous line : 
PNLF for the PNs bound to the M87 halo. Short dashed black line: Iciardullo et aD Jl989[ l 
analytical formula for the PNLF with m* = 26.0. Right panel : LOS velocity distribution 
for the 14 PNs bound to the M87 halo. The systemic velocity of M87 is also indicated. 



